function [G,lambda,vs,vp,vr,a,b] = calc_phy_params(rho,E,nu)
%calc_phy_params Summary of this function goes here
%   Detailed explanation goes here
G = E/(2*(1+nu));   % shear modulus
lambda = 2*nu*G/(1-2*nu);  % lame constant                                                                                                                                      

vs = sqrt(G/rho);   % s-wave velocity
kps = sqrt((2-2*nu)/(1-2*nu));          % kps = vp/vs
krs = -0.07602*nu^2+0.2012*nu+0.8738;	% krs = vr/vs
krp = krs/kps;      % krp = vr/vp
vp = vs*kps;        % p-wave velocity
vr = vs*krs;        % Rayleigh wave velocity
%
a = sqrt(1-krp^2);
b = sqrt(1-krs^2);
end

